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We formulate the conditions under which the dynamics of a continuously measured quantum 
system becomes indistinguishable from that of the corresponding classical system. In particular, 
we demonstrate that even in a classically chaotic system the quantum state vector conditioned 
by the measurement remains localized and, under these conditions, follows a trajectory charac- 
terized by the classical Lyapunov exponent. 
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The emergence of classical chaos from quantum me- 
chanics is probably the most important theoretical prob- 
lem in the study of the quantum to classical transition. 
Because of the absence of chaos in isolated quantum sys- 
tems W and the noncommutativity of the twin limits 
h ^ Q (the semiclassical limit) and t ^ oo (the late- 
time limit, necessary to describe chaos), the fundamen- 
tal mechanism of how classical chaos arises from quan- 
tum mechanics remains to be elucidated. While there 
has been much progress recently in the development of 
sophisticated semiclassical methods for chaotic dynami- 
cal systems [0, attempts to unambiguously characterize 
notions of chaos in the exact quantum dynamics ||3[ and 
to extract classical chaos as a formal semiclassical limit 
have been less successful: a rigorous quantifier of 'quan- 
tum chaos' on par with the classical Lyapunov exponents 
has yet to be found. And, since formal techniques have so 
far not succeeded in extracting trajectories from isolated 
quantum systems, they have not been able to explain the 
generation of chaotic time series in actual experimental 
situations. The experimental state of the art has, how- 
ever, reached the stage where the quantum to classical 
transition can now be probed directly |j]. So, it is cru- 
cial that one understand the mechanism underlying this 
transition in order to interpret existing results and design 
future experiments. 

As real experiments always deal with open systems, 
and the interaction with the measuring apparatus neces- 
sary to deduce classical behavior provides an irreducible 
disturbance on the free evolution of the quantum sys- 
tem, the resulting decoherence and conditioned evolution 
could play a crucial role in the emergence of the classical 
limit, and of chaos, from the underlying quantum dynam- 
ics. Indeed, some qualitative results in this direction al- 
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ready exist p, H| . In this Letter, we show that, even in the 
absence of any other interaction with the environment, 
the theory of continuous quantum measurements applied 
to the quantum dynamics of classically chaotic systems 
provides a quantitatively satisfactory explanation of how 
classical chaos, and Lyapunov exponents characterizing 
it, emerges from quantum mechanics. 

Open quantum systems are often studied by writing 
the evolution equation for the reduced density matrix ob- 
tained by tracing over the degrees of freedom in the envi- 
ronment. Although this process, called decoherence, can 
be extremely effective in suppressing interference effects 
and thereby making the quantum Wigner function ap- 
proach the corresponding classical phase space distribu- 
tion function p, it does not succeed in extracting local- 
ized 'trajectories' from the quantum dynamics. Without 
the existence of such trajectories it is extremely difficult, 
if not impossible, to rigorously quantify the existence of 
chaos both mathematically and in actual experimental 
practice. Since in order to extract classical trajectories 
systems must be observed, one expects observed quan- 
tum systems to obey classical dynamics in the macro- 
scopic limit. 

What is therefore desired is an unraveling of the Master 
equation which provides a more detailed understanding 
of the trajectories underlying the average system dynam- 
ics. When these detailed trajectories follow classical dy- 
namics (albeit noisy), one can infer that the average dis- 
tributions generated by them also become classical. This 
then provides a 'microscopic' understanding of the quan- 
tum to classical transition demonstrated, e.g., in Ref. [||. 

The first requirement in this program is to have a 
good model of continuous quantum measurement. Even 
though 'continuous' measurement is always an idealiza- 
tion, real experimental situations exist which approxi- 
mate it extremely closely, and simple models which corre- 
spond accurately to these processes have now been devel- 
oped 0, ^, |[ |I^ . These models show that as a necessary 
result of the information it provides, continuous mea- 
surement produces and maintains localization in phase 
space. On the other hand, the Ehrenfest theorem guaran- 
tees that well-localized quantum systems effectively obey 



classical mechanics. As the measurement process, in ad- 
dition to localizing the state, also introduces a noise in 
its evolution, to obtain classical mechanics one must be 
in a regime in which the localization is sufficiently strong 
and, yet, the resulting noise sufficiently weak. We show 
that such a regime exists, and is precisely the one which 
governs macroscopic objects, i.e. 7i <C S, the action of 
the system. In what follows, we refer to this regime as 
the classical regime. Our central result is that, once this 
regime is achieved, the localized trajectories for the con- 
tinuously observed quantum system obey the classical 
dynamics (possibly chaotic) for that system driven by a 
weak noise. As a result, even at a finite but non-zero 
value of h, the quantum 'trajectories' possess the same 
Lyapunov exponents as the corresponding classical sys- 
tem. As one goes deeper into the classical regime with 
?i — > 0, one can make the noise progressively smaller by 
optimizing the measurement, and, in the limit, the in- 
trinsic classical Lyapunov exponents are recovered. 

In order for this mechanism to satisfactorily explain 
the quantum-classical transition, the following conditions 
need to be satisfied: (1) localization as discussed above, 
(2) suppression of measurement noise, (3) the actual 
value of the measurement strength should become ir- 
relevant, and (4) the measurement record (i.e., the ac- 
tual results of the continuous measurement process), suit- 
ably band-limited, should follow the classical trajectory. 
These conditions are studied in more detail below. 

We consider, for simplicity, a single quantum degree of 
freedom, with position and momentum operators denoted 
by X and P, evolving under an unperturbed Hamiltonian 
P'^ /2m + V^(^). Quantum mechanics then dictates the 
familiar Heisenberg equations of motion for these opera- 
tors: X = P/m, and P = -dxV{X) = F{X). Except 
in the limit ?i — *■ 0, a continuous observation with finite 
measurement strength does not localize either the posi- 
tion or the momentum completely. Nevertheless, we can 
describe the state of the particle in terms of the central 
moments of X and P and, anticipating the limit, assign 
it to a point in phase space given by the mean values {X) 
and (P). 

The most natural measurement to use is a continuous 
measurement of position, not only because this is often 
what is observed with mechanical detectors, but also be- 
cause real schemes for the continuous measurement of 
position, considered in the field of quantum optics, may 
be described very simply ||ll[. In addition, a continuous 
measurement of position is an unraveling of the thermal 
Master equation in the high temperature limit, so that 
results demonstrated for this case also apply to decoher- 
ence due to a weakly coupled, high temperature thermal 
bath. We stress however, that we do not expect the par- 
ticular measurement model to effect the results signifi- 
cantly; any measurement or interaction which produces 
a localization in phase space should lead to classical be- 
havior in essentially the same manner. 

Under continuous position measurement the evolution 
of the wavefunction becomes stochastic. The stochas- 



tic master equation for the density matrix p(i), con- 
ditioned on the measurement record {X) + ^(<) with 
i{t) = {8Tjk)-^/'^dW/dt, is 



pit + dt) 
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2pTipX) dW, (1) 



where fc is a constant specifying the strength of the mea- 
surement, ry is the measurement efficiency and is a num- 
ber between and 1, and dW is a Weiner process, satisfy- 
ing (dW)"^ = dt. When ry = 1, the evolution preserves the 
purity of the state and can be rewritten in a way which 
allows it to be understood as a series of diffuse projection 
measurements WA on an unnormalized wavefunction ?/;: 
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where f (t), the difference between {X) and the measured 
value of the position, becomes a white noise in the limit 
that dt tends to zero. Under this continuous measure- 
ment process, the average values of position and momen- 
tum evolve according to 



d{X) ^{{P)/ni)dt- 
d{P) = {F{X))dt -+ 
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where V^ is the variance in position, and C^p is the sym- 
metrized covariance between X and P iQ. Thus the 
effect of the measurement is to provide some zero-mean 
noise proportional to the square root of the measurement 
strength k and to the width of the distribution. It is im- 
portant to note, however, that this is just the ffist in a 
hierarchy of equations for the moments; the equations 
governing the second moments contain terms depending 
on higher moments, and so on up the hierarchy. 

Even though the structure of the hierarchy makes it al- 
most impossible to obtain analytic answers to questions 
regarding the behavior of the variances, and resulting 
noise strength, which are the crucial quantities determin- 
ing the quantum to classical transition, we can, never- 
theless discuss the effect of varying h by truncating the 
hierarchy at the second order, and looking at the steady 
state solution for the variances. These equations show 
that to maintain enough localization to guarantee that, 
at a typical point on the trajectory, {F{x)) ~ F{{x)), as 
required in the classical limit, the measurement strength, 
fc, must stop the spread of the wavefunction at the un- 
stable pointsp^, dxF > 0: 
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Note that this condition is automatically satisfied for lin- 
ear systems, where quantum dynamics of the expectation 
values are identical with classical evolution. 

On the other hand, a large measurement strength in- 
troduces noise into the trajectory. If we demand that 
averaged over a characteristic time period of the system, 
the change in position and momentum due to the noise 
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FIG. 1: (a) The quantum trajectory in phase space, with 
h = 10~^ and k = 10^. (b) The position variance, Vx, as 
a function of time. 



are smah compared to those induced by the classical dy- 
namics, it is sufficient that, at a typical point on the 
trajectory, the measurement satisfy 



rjs 4 



(6) 



where s is the typical value of the action jlj] of the system 
in units of h. Obviously as s becomes much larger than 
2\/2?7~^'^, this relationship is satisfied for an ever larger 
range of fc, and this defines the classical limit. 

Finally, in experiments one usually considers the mea- 
surement record itself rather than the estimated state of 
the system as we have done. As measurement introduces 
a white noise, it is important to investigate the condition 
under which the record tracks the estimate faithfully. If 
At is the time over which the continuous measurement 
is averaged to obtain the record, and we allow ourselves 
a maximum of Ax as the position noise, it is easy to see 
that the measurement strength needs to satisfy 
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With this introduction, we consider, as an example, 
a bounded, one-dimensional, driven system with the 
Hamiltonian, 



H = P'^/2m + BX'^ - AX^ + AX cos(a)i) 



(8) 



with m = 1, B = 0.5, A = 10, A = 10, w = 6.07. This 
Hamiltonian has been used before in studies of quan- 
tum chaos [J_5| and quantum decoherence |5[| and, in the 
parameter regime used here, a substantial area of the ac- 
cessible phase space is stochastic. The numerical method 
used to solve Eq. (fi^) is a split-operator, spectral algo- 
rithm implemented on a parallel supercomputer. 

Simulations at various values of h confirm that as h is 
reduced, both the steady-state variance, and the resulting 
noise (for optimal measurement strengths) are reduced, 
as expected. As the dynamical time scale of this problem 
is 1 — 0.1, we decide to average the continuous observation 



record over a period of 0.01. Similarly, as the range of 
the motion covers distances of 0(10), we demand that the 
position be tracked to an accuracy of 0.01. By Eq. |j, this 
means we need rjk r^ O(IO^) or larger. In our example, 
we choose the energy to be O(IO^), and the correspond- 
ing typical action turns out to be 0(10), and the typical 
nonlinear! ty makes the rhs of Eq. g 0(1). We see that a 
choice of ?i = 10~^, rj = 1 and k — 10^, satisfies all the 
constraints for a classical motion. In Fig. H we demon- 
strate that in this regime, localization is maintained in 
spite of low noise. Fig. Il|(a) shows a typical phase space 
trajectory, with the position variance during the evolu- 
tion, Vx = (AA:)2, plotted in Fig. |I|(b). We find that the 
width AX is always bounded by 3.4 x 10^''. Further- 
more, as is immediately evident from the smoothness of 
the trajectory in Fig. 0(a), the noise is also negligible on 
these scales. 
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FIG. 2: (a) The quantum stroboscopic map with h = 
10^^ and k = 10^. The figure is a pastiche from several 
different runs with different initial conditions, for a total 
duration of 39, 000 periods of the temporal drive, (b) The 
stroboscopic map for the corresponding classical system, 
driven with a small amount of noise. 

Thus, we find that continuous measurement can effec- 
tively obtain classical mechanics from quantum mechan- 
ics. We substantiate this further by demonstrating that 
the trajectories we obtain show the common signatures 
of classical chaos. A direct way to compare qualitatively 
the global nature of the quantum and classical trajecto- 
ries in phase space is to compare the stroboscopic maps 
(the distribution of the locations of the system at a con- 
stant phase of the driving term). Fig. demonstrates 




FIG. 3: Lyapunov exponents, A, calculated for (a) the 
classical system driven with a small amount of noise, and 
(b) the continuously observed quantum system, with h = 
10~^ and k = 10^. The slope of the line drawn through 
the 
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cas e ' s is A ~ 0.57(2). 



the excellent correspondence between the classical and 
quantum maps in this regard. 

On a more quantitative level, we now calculate the 
key characteristic of chaos, the maximal Lyapunov ex- 
ponent, A, and compare it against that of the classical 
system driven with a similar amount of noise. We start 
with the definition of A: that for a chaotic system the 
distance between two nearby trajectories, A(i), evolves, 
on average, as InA(i) ^ Xt, as long as A(t) is small 
and t is large. To calculate this we take 10 fiducial tra- 
jectories (9 in the quantum case) starting at the point 
(-3,8) and at 17 points along each trajectory, separated 
by time intervals of 20 each, we obtain neighboring tra- 
jectories by varying the noise realization. The distance 
between the fiducial and these neighboring trajectories is 
tracked for a time interval of 8. The values of In A{t) thus 
obtained are averaged over all the instances and plotted 
versus t in Fig. pi both for the classical (with a small 
amount of noise) and the continuously measured quan- 
tum systems. For very small separations A, the noise 
dominates, which gives rise to an initial steep slope. This 
is followed by a linear region dominated by the Lyapunov 
exponent. Eventually A(t) becomes large and the curve 
flattens out. The behavior of the observed quantum and 
noisy classical systems are essentially indistinguishable, 
and the Lyapunov exponent 0.57(2) is the same for both. 
Performing the analysis with the classical system without 
noise, this time using 50 fiducial trajectories with initial 
points in a neighborhood of (-3,8), we obtain a Lyapunov 
exponent of 0.56(1), in agreement with the previous val- 
ues. 

After having demonstrated that in the classical regime, 
the localization and low noise conditions are satisfied si- 
multaneously, we study the sensitivity to the measure- 
ment strength. To this effect, we vary k between 2 x lO'' 
and 5 x 10^. The Lyapunov exponents remained un- 
changed within the quoted errors; only at fc = 5 x 10'"' did 



the noise start to wash out the flat region of the curve. 
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